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Abstract 

We consider neutralino production and decay e + e~ — > XiXji X°j ~^ Xiff 
at a linear collider with transverse e + and e~ beam polarization. We propose 
CP asymmetries by means of the azimuthal distribution of the produced neu- 
tralinos and of that of the final leptons, while taking also into account the 
subsequent decays of the neutralinos. We include the complete spin correla- 
tions between production and decay. Our framework is the Minimal Super- 
symmetric Standard Model with complex parameters. In a numerical study 
we show that there are good prospects to observe these CP asymmetries at 
the International Linear Collider and estimate the accuracy expected for the 
determination of the phases in the neutralino sector. 



1 Introduction 



Supersymmetry (SUSY) is at present one of the most prominent extensions of the 
Standard Model (SM) £Q. It can solve the hierachy problem, allows the unification 
of the gauge couplings and has the additional merit of providing new sources of CP 
violation. In the chargino and neutralino sectors of the Minimal Supersymmetric 
Standard Model (MSSM), the higgsino mass parameter \i and the gaugino mass 
parameter M\ are in general complex, while the SU (2) gaugino mass parameter 
M 2 can be chosen real by redefining the fields. The precise determination of the 
underlying SUSY parameters will be one of the main goals of the high-luminosity 
e + e~ International Linear Collider (ILC) 

The phases of the complex parameters \x and Mi may be constrained or correlated 
by the experimental upper bounds on the electric dipole moments (EDM) of electron, 
neutron and the atoms 199 Hg and 205 T1 These constraints, however, are rather 
model dependent. In a constrained MSSM the restrictions on the phases of /i and 
Mi can be rather severe. However, there may be cancellations between the different 
SUSY contributions to the EDMs, which allow larger values for the phases (for 
reviews see, e.g. [3]). If /i and Mi are complex and all other parameters are real, 
in general the phase of \i has to be small. However, the restrictions on the phase 
of \x may disappear if also lepton flavour violating terms in the MSSM lagrangian 
are included [2]. It is, therefore, necessary to determine in an independent way the 
phases of the complex SUSY parameters by measurements of suitable CP-sensitive 
observables. Experiments with transverse e beam polarization may allow us to 
construct suitable observables for precision studies of the effects of new physics and 
CP violation. Recent studies on the advantages of transversely polarized e + and e~ 
beams are presented in [13 HE]. 

The study of neutralino production 

e + e-^x-X?, i,j = l,...,A, (1) 
and subsequent two-body decay processes 

$ - - f^X? , n = l,2, (2) 

with £ = e, /i, will play an important role at the ILC. The production process has 
been studied extensively in the literature (see [2] and references therein). Production 
and subsequent decay processes, and the decay angular and energy distributions have 
been studied in detail in [H] J3]- The properties of Majorana and Dirac particles 
and their production and decay amplitudes under CP and CPT have been studied 
in [I2|-[in|- 

In [IB] JE] h has been shown that the parameters of the chargino and neutralino 
systems can be determined by measuring suitable CP-even observables. However, 
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the measurement of CP-odd observables is necessary to clearly demonstrate that CP 
is violated and for the unambiguous determination of the CP-violating phases. Sev- 
eral T-odd observables in neutralino production and decay applying triple product 
correlations have been proposed in [THl ED] • 

Transversely polarized beams offer the possibility to construct further CP-sensitive 
observables. This is the subject of the present paper. We propose and study CP-odd 
asymmetries for the case of transverse e + and e~ beam polarizations. We define two 
types of CP asymmetries, one that involves only neutralino production, the other 
one involving production and decay if one of the neutralinos decays as in Eq. (J2J). 
Note that in the case of chargino production it has been shown that the analogous 
CP-odd asymmetries vanish [21]. In neutralino production it is possible, however, to 
construct CP-odd asymmetries involving the transverse beam polarization, because 
there are t-channel and u-channel contributions due to the Majorana nature of the 
neutralinos jjj. The formulae for the production cross section of process (JTJ), for 
longitudinally and transversely polarized beams, have been given in In this 

paper we derive the compact analytic formulae for the CP asymmetries with the 
help of the spin density matrix formalism , including the complete spin correla- 
tions between production and decay of the neutralinos. We study numerically the 
parameter and phase dependences of the CP asymmetries. 

The paper is organized as follows: in Section |21 we set up the definitions. We 
present the formulae for the cross section of JTJ with transverse beam polarization 
in Sectional In Section 0] we define the CP-odd asymmetries. We present a numer- 
ical investigation of these asymmetries in Section while Section El contains our 
conclusions. 



2 Lagrangian and couplings 

The tree-level Feynman diagrams for the production process are given in Fig. ^ 
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Figure 1: Feynman diagrams of the production process e + e 
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The interaction Lagrangians are P 



C 



Z°£+£~ 



9 



l [L t P L + RtP R ]£ 



1 



cos Qw 



9 ry ~0..uro"£ d i r^'iRn l,~,0 



with z,j 



2 cos 0jy 

■ gff i lP R x% + gf i «£P L x% + h.c. 

4 and the couplings 



1 



cos 



11 



- + sin 2 Q w 



fg = V2sme w [t a ne w N* 2 - TV*] , 
O'lf = ~(N i3 N* 3 - N t ,N* A ) cos 2{3 - \{N iZ N* A + N^N* 3 ) sin 2/? 



O 



L e = -- + sin Q w , 



R f = sin 2 6 



(3) 

(4) 
(5) 

(6) 

(7) 
(8) 
(9) 
(10) 



where Pl,r = §(1-F7s)- is the weak coupling constant (g = e/ cos0vi/, e > 0), 0^ 
is the weak mixing angle and tan/3 = i^/fi is the ratio of the vacuum expectation 
values of the Higgs fields. The unitary (4 x 4) matrix iV^ diagonalizes the complex 
symmetric neutralino mass matrix in the basis (7, Z, H®, H®) |9 . 



3 Cross section 

In order to calculate the squared amplitude of process (JTJ) with subsequent decay 
(j2*j). we use the spin density matrix formalism j2H|. The squared amplitude of the 
combined process of production and decay reads 

|T| 2 = £ ^ |A(^)| 2 |A(x;)| 2 Pp iAj ' A:A; p d , a;Ai P aa;Aj , (11) 

A* At A - A ■ 

' J I J 

where A(%? •) is the propagator of the corresponding neutralino; pp 3 ' ! 3 denotes 
the spin density matrix of the production, p D x < A are the spin density matrices of 

the decay; Xij denotes the helicity of the neutralino Xij- The propagators are given 
by 

A(x° fc ) = l/[p|, - m\ k + im^ k ] , k = i,j. (12) 
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Here p^ k , m^ k and T^ k denote the four- momentum, mass and total width of the 
neutralino xt- For these propagators we use the narrow- width approximation. The 
(unnormalized) spin density production matrix is given by 

ii e e 

A e -A e+ A e _A e+ 

where Tp X 3 _ x is the helicity amplitude of the production process and A e ± is the 

' e e~i~ 

helicity of e ± . The spin density decay matrices can be written as 

Pd^x, = T dxT*d,x' z , ( 14 ) 
Pd.x'.x, = T DtX T* DX , , (15) 

3 J J 

where Td ^. is the helicity amplitude for the decay. The spin density matrices of 
the polarized e + and e~ beam in Eq. (|13|) can be written as 

P(e ± ) = \(l + Vta>) , (16) 

where Vf is the degree of transverse e ± beam polarization in the production plane, 
T>2 is the degree of transverse e ± beam polarization perpendicular to the production 
plane, = Vf; (—1 < Vf; < 1) is the degree of longitudinal beam polarization 
and a 1 (i = 1,2,3) are the Pauli matrices. For the degrees of transverse beam 
polarizations we have the relation (Vt) 2 + (Pt) 2 — (Pt) 2 with V 1 = cos (f)±P^ and 
Vt = sin0±P^ (0 < V$ < 1 and (p£) 2 + (P£) 2 < 1), see Fig. H 1 . 

The production and decay matrices are calculated with the help of the Bouchiat- 
Michel formulae 

u(p,X')u(p,X) = i[5 AA , +75 /V^](^ + m) , (17) 

v(p,X')v(p,X) = I[5 A , A + 75 ^J(^-m) . (18) 

The three four-component spin basis vectors s a and the 4- vector p/m form an or- 
thonormal system. For the incoming particles e + and e~ in the limit of vanishing 
electron mass, m e — > 0, Eqs. (|T7|) and (fTHjl can be written as 

u(p e -,X e -)u(p e -,X' e -) = -{(l + 2A e -7 5 )*V A +75[fe-4 A +fr4 A 

(19) 

: At this point we note that, contrary to the usual conditions at a circular accelerator, where the 
Sokolov-Tcrnov effect orientates automatically both transverse polarization vectors cither parallel 
or antiparallel (depending on the sign of the charge of the incoming particle), there is the possibility, 
at the ILC, to choose an arbitrary transverse polarization for both e + and e - , independent from 
each other. 
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V(j>e+, \' e +)v(p e +, A e +) = -{(1 - 2A e +7 5 )(V A +75[^ e 1 +^, + ^ e +C7^ ]} /f e + , 

Z e+ e e + e e+ 

(20) 

where t\± and tg± are the basis 4- vectors of the transverse polarization of the electron 
and positron beam, respectively. Thus, the transverse polarization 4-vectors can be 
written as 

t± = cosO± - <j>)t\± + sin(0 ± - <p)t 2 e± , (21) 

where is the azimuthal angle of the scattering plane and <p± are the azimuthal 
angles of the transverse polarization with respect to a fixed reference system, see 
Fig. |21 A convenient choice of t\± and t 2 e± in the laboratory system is given in 
Appendix |XJ 

With Eqs. (|17 j) - (j20| L the spin density production matrix and the spin density 
decay matrices can be expanded in terms of the Pauli matrices a a and a b , where the 
superscripts a (b) = 1, 2, 3 refer to the polarization vectors of Xi (x?) : 



Pp 3 = ^A^(X?«) + ^E^^(X?) 

a=l 

+Kx> E <.a' S p(x?) + E ^< A '.E?CX?X5) , (22) 

6=1 J a,b=l 3 

Pd,x% = W>(#) + IX* S Mx?) , (23) 



a=l 
3 



Paa'a, = ^^(x;) + E4,,^©- (24) 



The contribution P(x°x°) is independent of the neutralino polarization, whereas 
Sp(x°) and Ep(x°) depend on the polarization of the corresponding neutralino. 
Then Ep(x°j)/P(x°%°) gives the longitudinal polarization of the neutralino xljl 
Sp(x°j)/P(x°x°) is the transverse polarization of the neutralino in the scatter- 
ing plane; and Ep(x^)/P(x°x°) is the polarization perpendicular to the scattering 
plane. The terms Ep 6 (XiXj) describe the spin correlations between the polarizations 
of the two produced particles. The contribution of the decay matrix, which is inde- 
pendent of the neutralino polarization, is denoted by D(Xij), an d (Xij) denotes 
the contribution that depends on the neutralino polarization. 

With Eqs. [)22 |) -(|2% J) the squared amplitude |T| 2 , Eq. (fTTj). of the combined pro- 
cess of production and decay, for arbitrarily polarized beams, is given by 



|T| 2 = 4|A(x°)| 2 |A(x°)| 2 [P(^ 

CI=1 

+ ESp(x J )S 6 D (x J ) J D(x°) + E ^(x?X?)S!,(x?)ES,eg)] , (25) 



6=1 a,5=l 
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t X 




Figure 2: Decomposition of the polarization vectors into the longitudinal 
components Pl in the direction of the electron/positron momentum and the trans- 
verse components = V^t± with respect to a fixed coordinate system (x,y,z). 
The z-axis is in the direction of the electron momentum. 



where the contributions P(x°X°), £p(x°), Ep(x°), £p (x°x°) of the production spin 
density matrix include the terms for arbitrary beam polarizations. The differential 
cross section of the production and decay process is given by 

dcr = -^|T| 2 dLips(s,pi) , (26) 
2s 

where dLips(s,pi) = (2ir) 4 5 4 (pi + p 2 — J2iPi) Hi ^fy^E- ' see Appendix IE1 for more 
details. 



3.1 Contributions to the production spin density matrix 
independent of neutralino polarization 

In the following we regard the contribution -P(x°x°) to the spin density production 
matrix, see Eq. (}2*2"j) . that is independent of the polarization of the neutralinos x° 
and Xj- We only list the contribution PiXiX^T, which depends on the transverse 
e ± beam polarization: 

P(XiX j)T = P(ZZ) T + P{Z~e L ) T + P{Z~e R ) T + P(e L e R ) T , (27) 
P(ZZ) T = -V T V+ -J^-\MZ)\ 2 L e R e (\0'\ 2 + \0^)n , (28) 



P(Z~e L ) T = - V ^\^--R e 

2 cos^ <c>w 
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x {Re(A(Z)ftf i L J 0^[A(e L ,ty + A(e L ,uy}) n 

+ lm(A(Z)f e ff i L j O^[A(e L ,tr-A(e L ,ur}) r 2 } , (29) 



P(Ze R ) T = -VtV+ \ 



2 cos 2 ®w e 
x {Re(A(Z)ftf l pS[A(e R ,ty + A(e R ,u)*]) n 

- hn{A{Z)f?fjjOj}[A{e R ,ty - A(e*,ti)*]) r 2 } , (30) 



x {Re([A(e L ,t)A(e K , W )* + A(e L , U )A(e R ,0*]/^/^7#) n 

+ lm([A(e L ,t)A(e R ,uy - A(e L ,u)A(e R ,ty]f^ f^f f*) r 2 } , (31) 

where we have introduced the notation 

n = [(t-P4)(t+p 3 ) + (t-P3)(t+P4)}(piP 2 ) 

+[(pm)(P2P3) + {pm){p2Pi) - {piP2){pm)]{t~t+) , (32) 

r 2 = e^ pa [t + ^pi u p2pP4a(t-P 3 ) + t_^p lu p 2p P3a(t + p 4 ) 

+t- tl t +jU p 2p P4a{PlP3) + t-^t +>u p lp p 3a (p 2 P4)] , (33) 

where 5 0123 = 1 and A{Z) = i/{s - m|), A(e LjR ,t) = i/(t - mj LR ), A(e LjR ,u) = 
i/(u — "rrv~ L R ) with s = {j>\ + p 2 ) 2 , t = (pi — p^) 2 , u = (p± — pz) 2 - The masses of the 
selectrons and the Z-boson are given by rrie L R and mz- Since we study the process 
far beyond the Z-threshold, the Z-width can be neglected, and all propagators can 
be taken as purely imaginary. 

Note that only terms bilinearly dependent on transverse beam polarizations ap- 
pear for m e — > 0, because the couplings to e + e~ are vector- or axial-vector-like 
(23 I2H] (for the tL,R exchange the coupling to e + e~ can be brought to that form 
via Fierz identities [27]). Inspecting Eqs. (f2~%|) - ([3"Tj) . we note that transverse beam 
polarization gives rise to the interference term P(eLe R )T, which is absent for lon- 
gitudinal beam polarization [TI]. On the other hand, there are no terms P(eL<2L)T 
and P{& r &r)t for transversely polarized beams, but only for longitudinally polar- 
ized beams. Both are consequences of the Dirac algebra, since transverse beam 
polarization is described with an additional 7 matrix; see Eqs. (|19|) and (|2(J|) . 

The differential cross section for the process e + e~ —>■ is given by 

da=^^P(x?X?)dcos£d0, (34) 

where P (XiXj) contains the terms for arbitrary beam polarization and q is the 
momentum of the neutralinos in the center-of-mass system (cms) (see Appendix lA"|) . 
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3.1.1 CP-behaviour of the kinematical quantities 

In the following, the CP properties of the kinematical quantity r 2 , Eq. (|HHJ), and 
of the propagator difference [A(e, t) — A(e, u)] are discussed. These quantities con- 
tribute to the interference terms of the matrix element squared and are proportional 
to the imaginary parts of products of couplings, see Eqs. (f29|) - (j3l"j) . In the cms, r 2 
is given by: 

r 2 = 2E b [t + (fi x Pi)(t-p 3 ) + t-ipi x p 3 )(t + p 4 )} , (35) 

where Eb is the beam energy. Note that the second line of Eq. ()33|) vanishes in the 
cms. Applying a CP transformation to r 2 , Eq. (j3*3j) . with the following transforma- 

— * — » CP —> — » 

tions (pi,P2,f>3,P4,t~) < — ► (pi, P21 — P*3) — Pi: t+), we find that r 2 is CP-even. Since 
under CP A(e,t) « A(e,u) the propagator differences in Eqs. (|29| ) - (j3*T|) are CP- 
odd, their products with the CP-even quantity r 2 are CP-odd. We emphasize that 
this is due to the Majorana nature of the neutralinos, which leads to the simultane- 
ous presence of the t- and w-channel contributions, that the terms in Eqs. (|29 jl - (j31j) . 
which involve the imaginary part of the couplings, are non-vanishing in general. 



3.2 Contributions to the production spin density matrix 
dependent on neutralino polarization 

We now consider the terms Sp(x°) of the production spin density matrix, which 
depend on the polarization 4- vector s b of the neutralino In the following we 
only list the terms £p(x°)r, that involve the transverse beam polarization (for the 
contributions independent of the beam polarization and the terms that depend on 
the longitudinal beam polarization, see [TT]): 

S p (x?)t = Z b P (Ze L )T + X b P (Ze R ) T + E b P (e L e R ) T , (36) 



%{Z~e L ) T = V T V+l—^—R e 

2 cos z <c>w 



x {He(A(Z)/^/§Qg[A(g L) «)* - A(S L) f)*]) r b 



-lm(A(Z)ftfhO^(e L ,uy + A(e L ,ty]) r b } , (37) 



%{Z~e R ) T = V^\^—L e 

I cos^ k3w 



x {Re(A(Z)ftf e p^[A(e R ,uy - A(S H) f)*]) r\ 

+ Im(A(Z)/ £ f fgO$[A(e R , u)* + A(e R , t)*\) r b 2 } , (3? 
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x {Re([A(e L ,t)A(e R ,u)* - A(i L , u )A(e R ,t)*]ftf^ftfS) A 

+ lm([A(e L ,t)A(e R ,uy + A{e L ,u)A{e R ,ty]fI?fZf?fg) A} (39) 

with the following notation 

A = m Xj {[{t„ S b ){t + p 3 ) + {t„ P3 ){t + s b )){ Pl p 2 ) 

+ [(pis b )(p2P 3 ) + (Pip 3 )(p2s b ) - ( Pl p 2 )(p 3 s b )}(t-t + )} , (40) 

A = e^ vpa m X] [t + ^p lu p 2p s b a {t^p 3 ) +t^^p lu p 2p p 3a (t + s b ) 

+t_ fJ ,t +jU p 2p s b T (p 1 p 3 ) +t_ >fJ ,t +jU p lp p 3(T (p 2 s b )] , (41) 

where the e^ 1 polarization vector t± is given by Eq. (|2~Tf . The polarization basis 
4-vectors s b of the neutralino x°j fulfil the orthogonality relations s b ■ s c = —5 bc and 
s b -p± = 0. The parametrization of the neutralino spin vectors is given in Appendix 1X1 
The terms T, p (x®)t, which depend on the polarization 4- vector s a of x®, are obtained 
by the substitutions s b — > —s a ,m Xj — > m Xil p 3 — * p^ in Eqs. (HUJ) and (jJTJ). Note 
that, like P(XiX^)T, the expressions Ytpix^T contain no contributions Sp(ei,ei / )r 
and E b P (e R e R ) T , but an interference term, E b P (e L e R ) T . Furthermore, owing to the 
Majorana character of the neutralinos, there is no contribution H b P (ZZ)T- This is 
contrary to the cases of unpolarized and longitudinally polarized beams [TT] . 



4 CP asymmetries with transverse beam 
polarization 

4.1 CP asymmetries in neutralino production 

In this section we construct CP asymmetries for the production process e + e~ — > 
with transverse e + and e~ beam polarizations. The corresponding cross section is 
given in Eq. (|34j) . Choosing the e~ beam direction along the z-axis in the reference 
system (see Appendix El and Fig. |2J), the kinematical quantities in Eqs. fl3*2|) and 
(J2HJ) can be rewritten as 

n = -2El q 2 sin 2 9 cosfa - 20) , (42) 

r 2 = 2E 2 q 2 sin 2 9 sinfa - 20) , (43) 

where rj = 0_ + + . The CP-sensitive terms (oc r 2 oc sin(?7 — 20)) can be extracted 
from the amplitude squared by an appropriate integration over the azimuthal angle 
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We define the resulting asymmetry as 
A CP (9) - 



N[ 


sin (77 — 


20) > O;0] 


— N[sin(r/ 


-20) <0;£] 


N[ 


sin (77 — 


20) > O;0] 


+ N[sin(?7 


-20) <O;0] 


1 


\-r 




/•¥+! 




dV 


a 


V 

. 2 






2 — 2 . 


d0df? c 



(44) 



which depends on the polar angle 9. The first line in Eq. (J44J) exhibits how the 
asymmetry is obtained in the experiment, where N[sin(r/ — 20) > (< 0) ] denotes 
the number of events with sin(?7 — 20) > (< 0). The second line in Eq. (|44|) 
shows how the asymmetry is calculated. We can infer from Eqs. (l29|) - (p?Tj) that 
Acp{0), Eq. (|44[). would be zero if integrated over the whole range of 9: because 
of the propagators, the contribution of the t-channel cancels that of the u-channel. 
Therefore, we divide the integration over 9 into two ranges in order to obtain the 
CP asymmetry 



A 



CP 



sm(r] - 20) > 0;cos# > 0] - N[sin(r/ - 20) > 0;cos# < 0] 
+N[sin(77 - 20) < 0;cos# < 0] - N[sin(r/ - 20) < 0;cos# > 0]}/N 



tt/2 



tt/2 



A CP {9) M, 



(45) 



where N tot denotes the total number of events. Note that for a measurement of 
the CP asymmetry in Eq. (|45|) the production plane has to be reconstructed. In 
Appendix O we propose how this can be done. 

Finally we remark that an azimuthal asymmetry, analogous to that studied for 
chargino production 
by 



A A 



can be defined also for neutralino production. It is given 

N[ cos(r/ - 20) > ] - N[ cos(r/ - 20) < ] 
N[ cos(r/ - 20) > ] + N[ cos(r/ - 20) < ] 



a 







r — + 2 


r — + R 


r — + R ~ 










4^2 




/ -+ 


4 2 _ 


/ 4 2 + 










on ,v 










J — + 2 


4^2 . 



d0 # 



(46) 



In this case the integration over the polar angle 9 is performed over the whole range. 
This choice of the ranges of the integrations has the effect of extracting the terms 
oc ?"i oc cos(?7 — 20), Eqs. (|28 p — (|31| ) . from the squared amplitude. Note however, 
that this observable is CP-even. 



4.2 CP asymmetries in neutralino production and decay 



The reconstruction of the neutralino momenta is not necessary if we include the 
subsequent decays x°j — > ^ ± ^i (where £ = £l,£r) and £ ± - 



fx°i, yielding to the 
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final state Xj ~* ^f^txi- The label of the leptons indicates whether they stem from 
the first or the second decay. The cross sections for the combined processes are 
given in Appendix El Eqs. (J69)) and (fTUj) . respectively. The CP-sensitive terms of 
the squared amplitudes depend on sin(?7 — 2<^J or sin(r/ — 2</>£ 2 ), where 4>e.\ an d 
4>£ 2 are the azimuthal angles of the final leptons if and if. As a first step, we 
integrate the differential cross section in Eq. (jo^|) over all angles except <p£ l (the 
angles are integrated over their whole range). Then the CP asymmetry obtained by 
the azimuthal distribution of t[ is given by 

N[sin(r/-20 £l ) > 0] - N[sin(r/ - 2<j> tl ) < 0] 



N[sin(77-20 4 ) > 0] + N[sin(r/ - 2<j> tl ) < 



1 



+ 



+ 



2tt- 



3tt j y 
2 2 



do"i 



(47) 



where G\ = a(e + e — > XiXj) x B(Xj ~^ ) an d the upper index of A 1 corresponds 
to the electric charge of the observed lepton ly. 

As a next step, we integrate the differential cross section in Eq. ()70|) over all angles 
except 4>£ 2 , in order to define the CP asymmetry of the azimuthal distribution of if: 

N[sin(r/-20 £2 ) > 0] - N[sin(r/ - 20 £a ) < 0] 



Af 



N[sin(77-2<^ 2 ) > 0] +N[sin(77 



1 

^2 



/ 2 T 2 


r+3 _ 
















. J 2 


"2 + 2 





<0] 

dcr 2 

d0 £2 



(4? 



where 02 



cr(e + e 



n^r) x s(£ 



+ 



Note that, since 

S£)(x°) for the two C-conjugate decay modes of x°j ~^ i^i^ differs only by a sign 
(see Eqs. (J65|) and (|67j)) the asymmetries with upper indices + and — are related 
by Af = — A^ , 2 = 1,2. In order to measure both asymmetries, Eqs. (|47p and ()48|). 
it is necessary to distinguish the lepton originating from the decay x° — > i^if, 
and the lepton if from the subsequent decay — > Xi^f- This can be accomplished 
by their different energy distributions, when the masses of the particles involved are 
known, provided that their measured energies do not lie in the overlapping region 
of their energy distributions 



However, we can also define a CP asymmetry where it is not necessary to distin- 
guish whether the leptons stem from the first or the second step of the decay chain 
ififxi- This asymmetry is defined by 



0± 



A' 



N[sin(?7 



p-) > 



N[sin(7? - 2<f> e -) < 



N[sin(r/ - 2<f> t -) > 0] + N[sin(r/ 

a + -r)(^#4 



< 0' 



da2 , 

d<f>e 2 



(49) 
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where l~ is either £j or l 2 an d N[sin(?7 — 20^-) > (< 0) ] denotes the number of 
events where sin(?7 — 20£-) > (< 0). Hence, only the charge of the lepton and 
its azimuthal angle 4>e~ nas to be determined. In Eq. ()49)1 ^ corresponds to an 
integration over the azimuthal angles <pe 1 or <p£ 2 , where sin(?7 — 20^ 2 ) is positive or 
negative, respectively. An analogous asymmetry can be defined for l + as well. The 
asymmetry A~, Eq. (JUJ), can be related to the asymmetries A± and A 2 , Eqs. (}4T)) 
and (111, by 



A' 



1 



[1 + B{t 



-*?)] 



(50) 



5 Numerical studies 



5.1 CP-even observables in neutralino production 

Before we concentrate on the numerical study of CP-odd observables, we would like 
to give an example, which shows that a measurement of only CP-even observables 
may not be sufficient to unambiguously determine the SUSY parameters of the 
neutralino sector. However, a measurement of a CP-odd asymmetry may help to 
single out the correct solution. This may be particularly important if only the two 
lower states of the neutralino spectrum are kinematically accessible. 



Scenario 


|Mi| 




M 2 






tan/3 


m~ eL 


me R 


Complex 


183 


0.057T 


311 


343 


1.97T 


2.1 


297 


181 


Real 


180 





310 


335 





3 


300 


180 



Table 1: Input parameters \Mi\, 4>Mi, M 2 , rrie L and rrie R for the complex and 
the real scenario. All mass parameters are given in GeV. 



To this end we consider the complex scenario with the parameters given in Ta- 
ble d leading to m^o = 170.9 GeV and m^o = 259.5 GeV. At ^ = 500 GeV 
only the cross sections of e + e~ — > X1X2 would be measurable, giving <r(e + e~ — > 
X1X2) = (16.4, 18.3, 30.3) fb for the e + and e~ longitudinal beam polarizations 
(Pl^l) = (0,0), (-80%, +60%), (+80%, -60%), respectively. We assume that 
the masses m^o 2 and measured with 1% accuracy. For the cross sections 

we take an error corresponding to a 1-a deviation for a luminosity £; nt = 100 fb -1 . 
Then within this accuracy, we would obtain compatible neutralino masses and cross 
sections with the real SUSY parameter set, which is also given in Tabled namely 
m^o = 169.3 GeV and m^o = 258.3 GeV, and cross sections a(e + e~ — > X1X2) = 
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(a) e+e" - xlxl 
0.1 p — ' — ■ — ■ 1 — 




-0.1 V, r 

1 2 



<t>M 1 [ 7T ] 

Figure 3: (a) CP asymmetry A C p, Eq. 
as a function of 4>Mi for scenario A of TabJ2l 
(dashed line), tan/3 = 30 (dotted line), for 
polarizations (Vt,V£) = (100%, 100%). 



(b) e+e- - x?X^ 

200 P ' — 



180 - 




100 U ......... r 

1 2 



4>M t [ 7T ] 

and (b) cross section a(e + e~ — > X1X2) 
with tan/3 = 3 (solid line), tan/3 = 10 
y/s = 500 GeV and transverse beam 



(16.3,18.2,30.0) fb. The CP-odd asymmetry Acp, Eq. however, would re- 

sult in about 2.8% with (V^,V^) = (80%, 60%) for the complex scenario with 
0Mi = 0.057T and 0^ = 1.97T. Although the asymmetry is small it should be exper- 
imentally measurable including the statistical uncertainty. Therefore the complex 
scenario would be clearly distinguishable from the real scenario, which results in an 
asymmetry identical to zero. This simple example illustrates that it is necessary to 
measure CP-odd observables for truly identifying CP-violating effects. 

In the following we analyse numerically the CP-odd asymmetries, Eq. (J45[) and 
Eqs. P7)) -P9 |) . at the ILC with y/s = 500 GeV and transversely polarized 
beams. We especially focus on the influence of the phase <f>Mi of the gaugino mass 
parameter M\ = \Mi\e lr ^ M ^. Throughout we assume the GUT- inspired releation 
I Mi I = 5/3 tan 2 Qy/ M 2 . Furthermore we show that CP-odd observables are nec- 
essary to determine unambiguously the underlying SUSY parameters. In order to 
study the full phase dependences of the CP-odd observables, we do not take into 
account the restrictions from the EDMs and vary </> M and 4>m 1 in the whole range. 



5.2 CP-odd asymmetries in neutralino production 

First we discuss the CP-odd asymmetry Acp, Eq. (|45j). for the neutralino production 
processes e + e~ — > X1X2 anc ^ e+e ~ ~ * XiX®- CP violation is due to the interference 
terms P^Ze^T, P(Zeji)T, and -P(eie^)r, Eqs. (|2l?]) -([3*T |) . We assume that the mo- 
menta of the produced neutralinos can be reconstructed by analysing the subsequent 
two-body decays; see Appendix O 
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Scenario 


A 


B 




Mi 




123.3 


120.8 


4>Ah 


0.57T 


0.57T 


M 2 


245 


240 








160 


300 










m iL 


400 


400 


m~ eR 


150 


150 


tan/3 


3 


30 


3 




99.4 


105.5 


117.0 


X9. 


143.0 


144.1 


197.6 




169.7 


178.6 


303.9 


m xl 


289.7 


281.5 


351.7 



Table 2: Input parameters \Mi\, M 2 , |/i|, rrie L and m,e R and the resulting masses 
m^o, % — 1, . . . , 4 for tan (3 = 3, 30 and specific values of the phases 0Mi and 0^. All 
mass parameters and masses are given in GeV. 



5.2.1 CP-odd asymmetries in e + e — > X1X2 production 

In Fig.EK we show Aqp, Eq. (|43j) . for e + e~ — > X1X2 as a function of <f>M 1 for scenario 
A, given in Tabled for tan/3 = 3, 10, 30, with ^ = 500 GeV and transverse beam 
polarization (Vtj'Pt) = (100%, 100%). For this scenario we obtain for tan/3 = 3 
(30) an asymmetry A C p of about 8.2 (7.8)%, for <f) Ml = 0.5tt. The peculiar shape of 
the curve is a result of the combined contributions of the Z-e.R interference term, 
which has its maximum at 0^ pa OAn, and of the e^-e^ interference term, with 
its maximum at 0^ ~ 0.87T. The contribution of the Z— ei interference term is 
suppressed because of the large mass of the left-handed selectron. The cross sections 
for the process e + e~ —> X1X2 are plotted in Fig. Eb and are about 163 (144) fb for 
0Mi = 0.57T. Note that the cross sections are independent of the transverse beam 
polarization, because these contributions depend on cos 20 (sin 20), see Eq. f!42|) . 
and disappear if integrated over the whole range of 0. In Figs. OK and b we can 
clearly see the antisymmetric dependence of the CP asymmetry and the symmetric 
behaviour of the cross section on the phase <Pm 1 . It is therefore obvious that both 
kinds of observables are needed for an unambiguous determination of the phase. 
Note that Aqp can be sizeable even for values of 0Mi close to and ir, which would 
be favoured by the EDM constraints. 

Now we estimate the observability of the asymmetry. One assumes that the same 
degree of transverse beam polarization is feasible as for the longitudinal polarization 
(Py = 80% and V? = 60%). Since the CP asymmetry A C p depends bilinearly on 
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(a) A CP [%} for e+e- - xlxl (b) Ant [fb^ 1 ] for e+e~ - x°x° 




100 200 300 400 500 100 200 300 400 500 



[GeV] [GeV] 

Figure 4: (a) contours of the CP asymmetry Acp, Eq. (jUJ), in % for the process 
e + e~ —> X1X2 i n the |/i|-M 2 plane. The MSSM parameters are 4>m 1 = 0.57T, M = 0, 
tan/3 = 3, = 400 GeV and mg R = 150 GeV at y/s = 500 GeV with transverse 
beam polarizations (V^,V^) = (100%, 100%). (b) shows the contours of the lumi- 
nosity £i nt , Eq. (joTj) . needed to measure the CP-odd asymmetry Aqp at the 5-a 
level with degrees of transverse polarization (V? ,V£) = (80%, 60%). The light-grey 
region is experimentally excluded by the exclusion bound m~± < 104 GeV |25] . 



the degrees of transverse beam polarization V? (V?) of the e~ (e + ), see Eqs. fl2T?j)- 
(jHIJ), we have to multiply the asymmetry Aqp for (Vt,Vt) = (100%, 100%) with a 
factor 0.48. The luminosity C- m t required for a measurement with specific significance 
can be estimated as 

C mt = (K) 2 /[A 2 CP a] , (51) 

where M a denotes the number of standard deviations and a the corresponding cross 
section for neutralino production. We obtain a luminosity C- mt ~ 99 (124) fb _1 
needed for a discovery with 5-cr, for tan/3 = 3 (30) and 0m x = 0.57T. 

Figure E^i shows the contour lines of the CP asymmetry Acp, Eq. (|45|). at = 
500 GeV for e + e~ — > X1X2 i R ^e |/i|-M 2 plane. The MSSM parameters are chosen 
to be <f>Mi — 0.57T, cj)^ = 0, tan/3 = 3, m e - t = 400 GeV and m e - fl = 150 GeV. The 
largest CP-odd asymmetry Acp is attained for sizeable gaugino-higgsino mixing. 
If the beams are fully transversely polarized, (V^,Vt) = (100%, 100%), then Acp 
could reach up to about 8.8% for M 2 w 240 GeV and ss 140 GeV. With a higher 
centre-of-mass energy y^s = 800 GeV, the asymmetry Aqp increases to about 12%, 
because the cross section, which is the denominator of Acp, decreases stronger 
than its numerator. In this region of the parameter space the Z-cr interference 
term, Eq. (|30[). is the main contribution to the asymmetry Aqp- In a gaugino-like 
scenario, for instance M 2 = 250 GeV and = 450 GeV, the ei-ep interference 
term is dominant and the others are suppressed. Generally the e^-ep contribution 
to the asymmetry is small for e + e~ — > X1X2 an d Aqp is therefore reduced to about 
1.6%. In order to obtain a larger ep-ep contribution, a larger mass splitting of ep 
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(a) A CP [%] for e+e- -+ xlxl (b) A CP [%] for e+e" -> xgxg 




0.5 1 1.5 2 0.5 1 1.5 2 



0*1 W ^ [ 7T ] 

Figure 5: Contours of the CP asymmetry Acp, Eq. (jUJ), in % for the process (a) 
e + e~ — > X1X2 and ( D ) e + e~ — > X1X3 m the M -</>Mi plane, for scenario A with 
tan/3 = 3, see Table El for yfs = 500 GeV and transverse beam polarizations 
(Pf,P$) = (100%, 100%). 

and e R is necessary. If Trip T ~ mg R the interference term P{6lCr)t is very small, 
see Eq. ()31|) . In Fig. HJd we plot the corresponding luminosity £ int , Eq. (joTj) . for 
transverse beam polarizations of {Pj-iPt) = (80%, 60%). For the maximum value 
of Acp, a luminosity £ int of about 81 ftT 1 would be needed for a discovery with 5-cr. 

In Fig. we show the contour lines of Acp, Eq. (|45jl . for e + e~ — > X1X2 i n 
the 4> fl -4>M 1 plane for scenario A (Table EJ) at y^i = 500 GeV and with transverse 
beam polarizations (Pt,Pt) — (100%, 100%). We obtain a maximum value of 
the CP-odd asymmetry Acp of about 8.9% for M « 1.67T and 0a/! ~ 0.47T. In 
this scenario the 0^ and the M dependence are of the same order of magnitude. 
The main contribution to the CP-odd asymmetry originates from the interference 
term P(Ze R ) T , Eq. (|30|) . i.e. the primarily involved coupling is f t? fezQ 12- The 
corresponding cross section for ^ = 1.67T and §m x — 0.47T is about 139 fb and 
the luminosity £i nt for a discovery with 5-cr is about 99 fb -1 for transverse beam 
polarizations (Pt,Pt) = (80%, 60%). Note also in this case the CP asymmetry 
A C p can be sizeable for values of 0Mi and M close to and tt. 

5.2.2 CP-odd asymmetries in e + e — > X1X3 production 

Figure shows the contour lines of the CP asymmetry Aqp, Eq. (J45|) . in the (j)^- 
0Mi plane for e + e~ — > XiXs- As shown in the case before, a large gaugino-higgsino 
mixing is necessary to obtain sizeable CP asymmetries. We investigate scenario 
A of Table El at = 500 GeV and (7>f,7#) = (100%, 100%). In this scenario 
the maximum value of Aqp is about 9.8% for M rs 0.17T and 0m x ~ 1-27T. Again 
the main CP-violating contribution is due to the Z-e R interference term. In this 
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<t>M 1 [ 7T ] (j>M 1 [ 7T ] 

Figure 6: (a) CP asymmetry A C p, Eq. (J23j), and (b) cross section a(e + e~ — > X1X3) as 
a function of </>Mi for scenario A of Table with tan/3 = 3 (solid line), tan/3 = 10 
(dashed line), tan/3 = 30 (dotted line), for ^Js = 500 GeV and transverse beam 
polarizations (Vt,V£) = (100%, 100%). 

example the largest asymmetries are obtained for small values of For </> M = 0.17T 
and 4>m 1 = 1-27T the cross section a(e + e~ — > X1X3) is 76 fb and the luminosity for a 
discovery with 5-cr is about 150 fb -1 . 

Figure shows the CP asymmetry A C p, Eq. for the process e + e~ — > X1X3 
as a function of 0^ f° r scenario A defined in Table El for tan/3 = 3,10,30 and 
t/s = 500 GeV. For tan/3 = 3 (30) the asymmetry Aqp reaches its maximum of 
about 9.6 (7.4)% at <pM 1 — 1-25 (1.55)7r. Here again the dominant contribution to 
Acp comes from the interference term P(Zep)T', see Eq. (130)1 . Note that the maximal 
CP- violating phase 0m x = f (mod"7r) does not necessarily lead to the highest value 
of the asymmetry. The reason for this is an interplay between the 0Mi dependence of 
the cross section, shown in Fig. \§]p, and that of the numerator of the asymmetry. In 
Fig. Eh the corresponding cross section a(e + e~ — > X1X3) is plotted. For the maximal 
asymmetry it is about 78 (91) fb. In order to measure the asymmetry Acp at 5-cr, 
the required luminosity is C mt ~ 150 (217) ftT 1 . 

5.3 Neutralino production and subsequent two-body decays 

Now we discuss neutralino production e + e~ — > XiX'j with the subsequent decays 
X°j -> £r?T and £% -> Xi^t- We stud Y the CP-odd asymmetries, Eqs. (ET7l)-(EE9l. 
which are defined by the azimuthal distribution of the final leptons l\ and £2- In 
this case CP-violation effects arise from the contributions of the spin correlations 
of the decaying neutralino, Eqs. ()37 j) - (}3*9~j) . which depend on the transverse beam 
polarization. We give numerical examples for e + e~ — > X1X2 and e+e ~ ~~ > XiXl- 
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Figure 7: (a) CP-odd asymmetries Af (solid), Af (dashed) and A + (dotted), 
Eqs. (jUJ-fjinj), for the process e + e~ — > X1X2 ~~ > XiXi^f^t an d (b) the corresponding 
cross section as a function of hi scenario B (Tabled with tan /3 — 3. The centre- 
of-mass energy is fixed at ^/s = 500 GeV and the transverse beam polarizations are 
(p-,p+) = (100%, 100%). 



5.3.1 CP-odd asymmetries in e + e — > X1X2 production and decay 

In Figure we show the CP asymmetries Af 2 and A + , Eqs. (|47j ) — f)49j ) . as a function 
of for scenario B defined in Table 121 The beam energy is y^s = 500 GeV 
with degrees of transverse beam polarizations (V^,V^) = (100%, 100%). We study 
neutralino production e + e~ — > X1X2! with the subsequent decays X2 ~^ ^R^f an d 
if^ — > x?^- For ^1" we obtain a maximal value of about 11.6% for = 0.457T. 
The asymmetry Af is reduced to 2.1% by the additional contribution to the phase 
space from the decay l\ — > Xi^f- Since the branching ratio B(£^ — > = 1) 

we obtain a CP-odd asymmetry A + « 6.9%, see Eq. f|50|) . for 0Afi = 0.457T. In this 
scenario the main contribution to the CP-odd asymmetries comes from the el-en 
term, Eq. (|39j) . With a smaller mass splitting between and the contribution 
to the asymmetry of £p(e£,ep)p becomes larger, but as the cross section is increasing, 
the combination of the two effects leads to a smaller asymmetry. The corresponding 
cross section u(e + e~ — > X1X2 ~^ XiXi^^f) 1S plotted in Fig. 03. For <$>m x = 0.457T, 
we obtain a cross section of about 46 fb. Thus, for (V^,V^) = (80%, 60%), the 
luminosity £i nt needed for a discovery with 5-cr of the asymmetry Af is about 
176 fb -1 . For a discovery with 5-cr of A + , C mt pa 517 fb -1 are needed. 

In Fig. |HK we show the contour lines of the CP-odd asymmetry Af, Eq. (|47jl . for 
e + e~ — > X1X2 ~~ > Xi^r^i i n the |/x|-M 2 plane. The other parameters are 4>m 1 = 0.5/T, 
M = 0, tan/? = 3, m gi = 400 GeV and m gfl = 150 GeV. The centre-of-mass 
energy is fixed at ^/s = 500 GeV with degrees of transverse beam polarizations 
{Pti'Pt) = (100%, 100%). In this figure we only consider the parameter regions 
where the decay channel X2 — > ^r^x ls kinematically accessible. The maximum 
value of the CP asymmetry Af pa 12.6% is obtained for a gaugino-like scenario 
with M2 = 200 GeV and \fj,\ = 280 GeV. For these parameters the neutralino 
masses are m^o = 97 GeV and m^o = 163.5 GeV, and therefore the branching ratio 
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A+ [%] for A+ [%] for 
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Figure 8: Contours of the CP-odd asymmetry Af , Eq. ()47)l . in % in the |/i|-M 2 
plane for the process (a) e + e~ — > X1X2 — * Xx^R^t an d (b) e + e~ — ► X1X3 ~~ *• Xi^r^i- 
The MSSM parameters are 0a/! = 0.57T, </> M = 0, tan/5 = 3, m gi = 400 GeV and 
mg R = 150 GeV. The centre-of-mass energy is fixed at \fs = 500 GeV and the 
transverse beam polarizations are {Vt^t) = (100%, 100%). The light-grey region 
is excluded by m-± < 104 GeV |2E|. 

B (X°2 ^1) = 1 and the cross section a{e + e~ -> x\xl Xi^f) = 54 fb. 
Thus for (Vf,Vx) = (80%, 60%) the required luminosity £ int for a discovery with 
5-cr, Eq. (|51jl . is about 128 fb -1 . For this parameter point the asymmetry A + , 
Eq. is about 6.7% and the luminosity £i nt ~ 456 fb -1 . For higgsino-like 

scenarios (M 2 > \fi\) the Ep(Zefj) T interference term, Eq. (f3*5Jl. gives the main 
contribution to the CP-odd asymmetry, which can be traced back to the structure of 
the corresponding coupling f a* f^O 12- On the other hand for gaugino-like scenarios 
> M 2 ) the contribution of the interference term Sp(eie^)T dominates, with the 
corresponding coupling ftifkfiifei- The sign change of the asymmetry in the 
M 2 plane is therefore due to a cancellation of the Z-Sr and the cl-gr contributions 
which have opposite signs. 



5.3.2 CP-odd asymmetries in e + e — > X1X3 production and decay 

In Fig. [HJd we show the contour lines of the CP asymmetry Af for the process e + e~ — » 
X1X3 ~~ > Xi^R^t m the |yu|-M 2 plane. We fix the MSSM parameters at 4>Mi = 0.57T, 
M = 0, tan/5 = 3, m~ eL = 400 GeV and m iR = 150 GeV with yfe = 500 GeV 
and (Vf,V$) = (100%, 100%). The maximal CP asymmetry Af is about 31% for 
M 2 = 300 GeV and \fj,\ — 160 GeV. For this parameter point we obtain neutralino 
masses of m^o = 115 GeV and m^o = 156 GeV, and therefore the branching ratio 

is again -B(x 2 ~~ * = 1- Hence the cross section a(e + e~ — > X1X3 ~~ ¥ Xi^r^T) = 

83 fb. For transverse beam polarizations (V^,Vf) = (80%, 60%) the luminosity £ int , 
Eq. I|51|) . for a discovery with 5-cr of is about 14 fb -1 . For these parameters the 
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Figure 9: CP-odd asymmetries Af (solid), A\ (dashed) and A + (dotted), Eqs. (jUJ)- 
(|49|) . for the process e + e~ — * X1X3 — > XiXi^T^t an d (b) the corresponding cross 
section as a function of (pM 1 in scenario A, see Table with tan/3 = 3. The 
centre-of-mass energy is \fs = 500 GeV and the transverse beam polarizations are 
(p-,p+) = (100%, 100%). 



CP-odd asymmetry A + is 17.3% and the necessary luminosity (for a discovery with 
5-<r) Ant ~ 43 fb -1 . In the case of e + e — > X1X3 the CP-violating contributions 
Y^pi^ZeRjT and Ep(ei / e^)r, Eqs. (jHEJ) and (jSHJ), enter with the same sign due to 
the corresponding couplings. In gaugino-like scenarios the contributions of both 
interference terms are suppressed, therefore the asymmetry decreases. 

In Fig. ink the CP-odd asymmetries Af^ and A + , Eqs. (I47jl - (|49jl . for the pro- 
cess e + e~ — > X1X3 ~~ > XiXi^i^t are plotted as a function of <f>Mi for scenario 
A, see Table |21 For a centre-of-mass energy y/s = 500 GeV and with transverse 
beam polarization (Vf,V^) = (100%, 100%) the maximum of the CP asymmetry 
Af(A + ) w 26.2 (14.3)% is obtained for 0^ rs 0.757T. In this scenario we have large 
mixing between the gaugino and the higgsino components, the main contribution to 
Af stems again from the Z-e^ interference term, which is about 21.4%. The cr-cl 
contribution is 4%, wheras the Z-ei contribution is suppressed by the large mass of 
the left selectron. Figured shows the cross section a(e + e~ — > X1X3 ~^ Xi&R^t)- F° r 
the maximum of the CP asymmetries Af(A + ), for §m x ~ 0.75-7T, the cross section 
is 78 fb. For (Vt^t) = (80%, 60%) the luminosity C- m t for a discovery with 5-cr is 
about 20 (68) fb" 1 . 



5.4 Determination of the SUSY parameters 

In the following we will give an example for the accuracy that can be expected 
in the determination of the MSSM parameters, focusing on the determination of 
the complex parameter M\ = \Mi\e l< ^ M ^ . In order to determine the parameters 
unambiguously, CP-even as well as CP-odd observables have to be included in the 
set of observables from which the underlying parameters are extracted. 
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(0,0) 


(-80%, +60%) 


(+80%, -60%) 


a(e + e- -> xlxl) 


47.27 fb 


87.13 fb 


52.80 fb 


a(e + e~ -> X2X2) 


11.59 fb 


33.12 fb 


1.186 fb 




9.83 fb 


5.68 fb 


23.42 fb 


(j(e+e- -> 


7.86 fb 


7.74 fb 


15.53 fb 



Table 3: Cross sections for different sets of longitudinal beam polarizations in 
scenario B with 0Mi = 0.57T and M = for = 500 GeV. 

Our set of observables contains the neutralino masses m^o, the cross sections 

e + e~ — > XiX'j f° r different choices of longitudinal beam polarizations {Vi,V^) = 
(0,0), (-80%, +60%), (+80%, -60%) and the CP-odd asymmetry A CP , Eq. (jfejl . 
We now take scenario B with <pM l = 0.5tt and 0^ = 0, see Table El as our ref- 
erence point of input parameters. We calculate the neutralino masses m^o, m^o, 
m^o and m^o, see Tabled The cross sections for e + e~ — > X1X2, e + e~ — > X2X2' 
e + e~ — > X1X3 an d e + e~ — > XiXa f° r \A — 500 GeV with different sets of lon- 
gitudinal beam polarizations are displayed in Table 01 The CP asymmetry Aqp 
for the process e + e~ —>■ X1X2 * s about +2% for transverse beam polarizations 
(Vt,Vt) = (80%, 60%). We regard these calculated values as real experimental 
data, where we assume errors of 1% for the masses. For the error of the cross 
sections of each polarization configuration and of the asymmetry we take a 1-a de- 
viation for a luminosity £ int = 100 fb" 1 . Our approach for the determination of 
the error of the parameters (in particular the error of Mi) is described as follows: 
we perform a random scan over the input parameters \M\\, 0Mi, M2, \ \i |, M and 
tan/3 around our reference point and select the points which pass the condition 
|£,meas _ £)cak| < |A0f eas |, w h ere Qmeas are V8b \ ues Q f the observables at our 

reference point, see Table 01 AO^^ is the corresponding error, and (9^ alc are the 
values of the calculated observables obtained through the random scan. 

In Fig. we show the SUSY parameter points compatible with our reference 
scenario in the Re(M 1 )-Im(M 1 ) plane. If we consider only the CP-even observables 
(cross sections and masses) then we obtain two regions in the parameter space com- 
patible with our reference scenario. This ambiguity can be resolved if one includes 
in addition the CP-odd observable Aqp, Eq. (|45|). For the error of Mi we obtain: 
Re(Mi) =0 + 5.9 GeV and Im(Mi) = 120.8 ± 1.3 GeV. 

6 Conclusion 

We have studied the processes e + e~ — > X1X2 an d e + e~ — > X1X3 with subsequent 
decays X23 ~~ > ^R^ an d £r — > Xi^ where I — e, fi, at a linear collider with transverse 
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Figure 10: SUSY parameter points in the Re(Mi)-Im(Mi) plane consistent with 
scenario B, if one assumes an uncertainty of 1% for the masses and l-er deviation 
(Ant = 100 fb -1 ) for the cross sections and the asymmetry. The parameters \Mi\, 
0Mi, M2, I n |, (pfj, and tan/? have been randomly scanned around the reference point 
B. The grey points are excluded if one takes into account the CP-odd observable 
A C Pi Eq. (03). 



e + and e~ beam polarizations. We have discussed different CP asymmetries, which 
are due to azimuthal distributions of the neutralinos or the final leptons. We have 
pointed out that these CP asymmetries are non-vanishing thanks to the Majorana 
character of the neutralinos. We have given the analytical expressions for the CP 
asymmetries and the cross sections in the spin density matrix formalism, including 
the complete spin correlations between production and decay of the neutralinos. At 
the ILC at y^s = 500 GeV and with degrees of transverse beam polarizations 
{Pti'Pt) = (80%, 60%), the CP asymmetries can reach up to about 15%. Also 
we have shown that these CP asymmetries can be observed in a broad range of 
the MSSM parameter space. Furthermore, we have discussed the unambiguous 
determination of the underlying SUSY parameters, which requires CP-even as well 
as CP-odd observables. 
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Appendices 



A Momentum and polarization vectors 

The basis vectors for transverse e~ beam polarization are 

£- = (5- XPe-Vl^e- ><Pe-| , (52) 
= (Pe- XP Xj )/\Pe- Xp Xj \ ■ (53) 

The basis vectors for transverse e + beam polarization are defined analogously. In a 
fixed coordinate system (x, y, z), with the z-axis pointing along the beam direction 
the basis vectors in the cms are given by 

t\± = (0, cos0, sin 0, 0) and t 2 e± = (0, — sin0, cos0, 0) . (54) 

The momentum 4-vectors of Xi an d X°j are 

Pxj,v = P*,i* = q(E Xj / q, cos (f> sin 9, sin <j) sin 9, cos 9) , 

Pxi,» = P3,n = q(E x J q, - cos 4> sin 9,- sin 4> sin 9, -cos 9) , (55) 

with 

s + ml —ml.. \i (s, ml., ml.) 

E = = v ^ (56) 

where X(a, b, c) = a 2 + b 2 + c 2 — 2(ab + ac + be). The three spin-basis vectors s b x .^ 
of are chosen to be 



si... = I 0, I = (0, — cos d>cos9, — sin 6 cos 9, sin 9) , 

3 V P2 X S 3 \J 



\ ( E \ E 
s Xj fi = g, —^-p X j I — ~^~{q/E xp cos</>sin6>, sin^sin^, COS6 1 ) , (57) 



where s~ , s^ 2 and si build a right-handed-system. The momentum 4- vector of the 
lepton in the decay x° - ► ^l,_r^ is given by 

Pi!,v= \pe 1 \(l,cos4> h sin9 ei ,sin(j) h sin9e 1 ,cos9 ei ) (58) 

with 

2 2 
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and 

cosi? = sin # sin 6^ cos(0 — 0^) + cos # cos 6^ : . (60) 
The momentum 4-vector of the lepton from the decay I — > Xi^2 is given by 

Pi 2 ,fi = \Pta\0-, cos 04 sin 9 i2 , sin <f) h sin 9 h , cos 9 i2 ) , (61) 

where 

2 2 
mf — mi 

\p £2 \ = ^ . (62) 

B Decay matrix and phase space of 2-body decay 



The spin density matrix of the decay x° — > ^ can De written as 

Pd,x>x 3 = fy^Cx?) + E ^a,^(x°) , (63) 



c=l 

where the expansion coefficient -D(Xj) is the P ar t that is independent of the polar- 
ization of the decaying neutralino Xp an d S'MXj') is the part that depends on the 
polarization of x°j- For the sake of simplicity we consider I = e, /i, where the mixing 
in the slepton sector can be neglected. Then we have, for x°j ~ > 

D(xj - ll ^) = y l^l 2 K - m V > ( 64 ) 
^(x? - § ^) = ^ 2 |^| 2 m Xj (^ " Pi) (65) 

r,o _^ t>T e± 



and for x° — > £^ £ 



^(x? - ll = ±9 2 \fS\ 2 m X] (s c ■ p e ) , (67) 

where m Xj (tti^ l ) is the mass of x® (^l,r)- The parametrizations of the momentum 
4-vector and the polarization 4-vector s~. ^ of the neutralino x° are given in 
Eqs. (JS3j) and (JSTj) in Appendix EI Finally, the matrix element squared for the 
two-body decay £ T — > Xi ^ i n t ne decay chain, Eq. 0, is 

- X? O = £ 2 |/a'*| 2 K fl - O ■ (68) 
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From (|25|) and (|63|L summing over the polarization of Xi\ whose decay is not con- 
sidered, the differential cross section for e + e~ — > xlx'j ~~ > Xa^l.r ^1 is: 

dai = -[PD+ £ P £^] | A(x°)| 2 dLip Sl . (69) 
s J 

Similarly one obtains the differential cross section for e + e~ — > — ¥ Xi^L,R @t ~> 
ZftZXiXi from (EH), (El and flU) and summing over the polarization of Xi '■ 

da 2 = -[PD + D 2 |A(x°)| 2 |A(£)| 2 dLips 2 , (70) 

s 

where P and Ep involve the terms for arbitrary beam polarization. For the calcula- 
tion of the cross section we use the narrow widths approximation (/ |A(x°)| 2 ds^ = 
— ^fr- and J \A(£)\ 2 dSf = — ^p-, where ds$. = p%. and dsp = p~). The Lorentz- 
invariant phase-space elements in Eqs. (J69)) and (|7Uj) for the decay chain x° — > 
- X?^2 T are 

dLips! = —dLips(s,p Xi ,p x .)ds m dUps(s m ,p i ,p h ) , (71) 
dLips 2 = -^-pdLips(s, P Xi ,P Xj )ds mx . dLips(s mx . , p h p ei )ds m AUps(s mp p Xl , pt 2 ) 
with the Lorentz invariant phase space elements 



(72) 



1 q 

dLips(s, p Xt , p Xj ) = ^ sin d# d<p , (73) 

1 I ~* I 

dLips(s m , P/ , p 4 ) = — — - P£l sin 4 d0 £l d<fo , (74) 

A{A7i) m % . - m- 

dLips(s m p p xi ,pe 2 ) = J" , 2 fi 2 - 2 sin^ 2 d^ 2 d<^ 2 . (75) 



2(27r) 2 m?-m 2 



C Reconstruction of the production plane 

As an example, we consider the process e + e~ — > X1X2 with the decays x% ~~ > ^1 ^ 
and £ ^ £2 x'i-, where we denote the neutralino from the decay by x'i (here again 
the labels of the leptons indicate their origin). We assume that the masses of all 
particles involved are known. 

We rotate to a coordinate system where the 3-momentum vector of £\ is along 
the z-axis and that of £2 is in the x—z plane. The unit vectors of the 3-momenta of 

£1, £ 2 , £ are 

Pi x = (0,0,1) , p? 2 = (sine, 0, cose) , Pi — (sin b cos A, sin b sin A, cos b). (76) 
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From the relation (pe 1 + pi) 2 

cos b = 



p 2 2 we obtain 



1 



2 \Pei\\Pe 



\2 I 



Wh\ 2 - \n 



(77) 



where \p^\ 2 = E~ —mj and Ej> = E X2 — Eg 1 . From a second relation {p^- 



P 



between momentum vectors we obtain 

1 



xi 



cos a 



2|p< 3 |IhI 



(2 



(7J 



4 



m~ and i? v ' 



Ee 



where Ixv 1 2 = i? 2 , 

l-^Xi I Xi Xi Xi X2 

is the missing energy. From spherical geometry with p« ■ p~ t 
following relation between the angles 



or E v i 

Xi 



cos A 



cos a — cos b cos c 
sine sin b 



-E - E X1 , and E 
cos a, we obtain the 



(79) 



Inserting ((77)) , (J7HJ) and (|75jl into p| in (|7B)l. this vector is determined up to a twofold 
ambiguity in the second component. In order to resolve this ambiguity, a reference 
vector is needed, which tells us in which hemisphere of the x-z plane the momentum 
vector pi is. For instance, this is possible in the process e + e~ — > X1X3 with the decays 
X3 — > £1 £, I — > £ 2 X2 anc i X2 ~ ^ Xi -^; where the 3-momentum of the Z boson is the 
reference vector. 
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